3 May 
LPTENS-96/23 
UT-Komaba-96/11 



CORRELATIONS OF NEARBY LEVELS INDUCED 
BY A RANDOM POTENTIAL 

E. Brezin") and S. Hikami'') 



Laboratoire de Physique Theorique, Ecole Normale Superieure 

24 rue Lhomond 75231, Paris Cedex 05, France^ 
Department of Pure and Applied Sciences, University of Tokyo 
Meguro-ku, Komaba, Tokyo 153, Japan 



Abstract 

We consider a Hamiltonian H which is the sum of a deterministic 
part Hq and of a random potential V . For finite N x N matrices, 
following a method introduced by Kazakov, we derive a representation 
of the correlation functions in terms of contour integrals over a finite 
number of variables. This allows one to analyse the level correlations, 
whereas the standard methods of random matrix theory, such as the 
method of orthogonal polynomials, are not available for such cases. At 
short distance we recover, for an arbitrary Hq, an oscillating behavior 
for the connected two-level correlation. 
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1 Introduction 



Let us first recall the results for the correlations between two eigenvalues 
for the simple unitary ensemble, in which the full Hamiltonian is treated as 
random. In the simplest Gaussian ensemble one considers N x N random 
Hermitian matrices H with probability distribution 

PiH) = lexp(-^M2) (1.1) 

Two kinds of universal correlations between eigenvalues are known to be 
present for such problems: a) a short- distance universal oscillatory behavior; 
b) a finite distance universality of smoothed correlations. 

Let us review these two properties. The density of eigenvalues and the 
two-level correlation function are defined as 

p(A) =< ^Tr5(A - if) > (1.2) 

and 

p(^)(A,/i) =< ^Tr6iX-H)j^TrS{fi-H) > (1.3) 

The correlation function, when A and fi are arbitrary, has a complicated, 
non- universal, oscillatory behavior. It simplifies and is independent of the 
probability distribution of H, when 

a) A — /i is small, N is large, and the scaling variable 

x = 7riV(A-/i)p(^(A + /i)) (1.4) 



1 ,- / , N / , N / , N / N sin^ X 



is held finite. Then one findsUl] 

Pf\\p) ^ ^5(A-p)p(A)-p(A)p(/i) ^2 

~ WX^ 1 sin^[7riV(A-p)p(^)] 

- ^5(A-/.)P(A)-^ (1.5) 

b) Away from this short-distance region, for arbitrary A and p, the cor- 
relations simplify only if one "smooths" the oscillations. This is what one 
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usually does, if one lets N go to infinity first in the resolvent, before returning 
to the real axis. The result, which is known to be universal, is 

{2)Q ^ = 1 1 - V) ^^ 

Pc K^.N 2Ar27r2(^_^)2[(^2_^2)(^2_^2)]l/2 l''^^ 

where a is an end point of the support. 

There are many equivalent derivations of the property b). They are based 
either on orthogonal polynomials 0], or on summing over planar diagrams 
[0, or solving an integral equation however the property a) is known 

only through the orthogonal polynomials approach [0]. For the generalization 
that we have in view in this article, in which the "unperturbed" part of the 
Hamiltonian is deterministic, if again for b) a diagrammatic approach still 
works 0, ^, 0, H, we are not aware of any method which would allow us 
to study whether a) still holds. To this effect we shall generalize a method, 
introduced by Kazakov0, to the study of correlation functions. It consists 
of introducing an external matrix source. It leads to a representation of the 
correlation function in terms of contour integrals over two variables for finite 
N. We have used already this representation in a previous paper devoted to 
random matrices made of complex blocks, i. e. the Laguerre ensemble |1I0|] , 
but there we have let the source go to zero at the end of the calculation. 
Keeping this source finite allows one to deal with an arbitrary deterministic 
Hq. We shall illustrate it here when the random potential V belongs to the 
simple Gaussian unitary ensemble. 

2 An external matrix source: deterministic 
plus random hamiltonian 

We consider a Hamiltonian H = Hq + V , where Hq is deterministic and V is 
a random N x N matrix. The Gaussian distribution P is given by 
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= }_^-fTriH^-2HoH+H^) ^2.1) 

We are thus simply dealing with a Gaussian unitary ensemble modified 
by a matrix source A = —Hq. Up to a factor the probability distribution for 
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H is thus 

P^(H) = -^exp{-—TiH^ - NTtAH) (2.2) 

Let us first show how one deals with the density of states p(A) . It is the 
Fourier transform of the average "evolution" operator 

UAit) =< Iw*^ > (2.3) 

and p(A) is 

We integrate first over the unitary matrix u which diagonalizes H, and with- 
out loss of generality we may assume that A is a diagonal matrix with eigen- 
values (ai, ■ ■ ■ , otv) . This is done by the well-known Itzykson-Zuber integral 
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/<.e.p(T.W,^^|;^ (.5, 

where A (A) is the Van der Monde determinant constructed with the eigen- 
values of A: 

N 

A{A) = n(a. - a,) (2.6) 

i<j 

We are then led to 

1 1 ^ r 

X exp^-JY,'■"-'^^J2''f^) (2.7) 
The normalization is 

f/A(0) = 1 (2.8) 
The integration over the rj may be done easily, if we note that 

/ rfri • • ■ driyA^ri, rAr)exp(-y ^ ^Y. ^'^i) 

= A(6i,---,6^)efS^? (2.9) 
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If we use this, with 

it 

hi = ai- j^^a,i (2.10) 

we obtain 

UA{t) = ^ E n ( ""~!:~^ )e-^-"°" (2.11) 

The sum over N terms may be replaced by a contour-integral in the complex 
u plane, 

UA{t) = --J^I[ (^i^^l^)e---^ (2.12) 

The contour of integration encloses all the eigenvalues a-y. If we let all the 
go to zero, we obtain 

f/o(t) = -ie-w / ^e--(l - ^ f (2.13) 

From this exact representation for finite N, it is immediate to recover all the 
well-known properties, the semi-circle law, or the more subtle edge behavior 
of the density of states. Since the result is similar to the Laguerre case that 
we have discussed in an earlier paper |T^, we shall not discuss it here. 

For the two-level correlation function, p^'^\\, jj) is obtained from the 
Fourier transform UA{ti,t2), 

P^'W^) = I I ^e-*^^-*^^t/^(ti,t2) (2.14) 



where UAiti, ^2) is 



(27r) 



UAihM) =< ^Tre^^^^lTre^*^^ > (2.15) 



The normalization conditions are 

UA{tiM)=UA{t2M) 

f/^(ti,0) = t/^(ti) 

Ua{^) = 1 (2.16) 

Dealing with f/4(ti, ^2) is also simple; we have exposed the technique in more 
details in the almost similar problem of the Laguerre ensemble|]TO|. After 
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performing the Itzykson-Zuber integral over the unitary group as in 
we obtain through the same procedure, 

[/f(ti,t2) = ^ E In^^^lr^re"''^^^^-'^""^^^^^^^^ (2-17) 

cn,Q!2 = l i=l ^ ' 

After integration over the rj, we obtain 

^ aT:^2 Ui<j{ai-aj) 

f2 ,2 

^g— Itiact]^— «t2aQ2~2]V^2JV JV^'''»l'<»2 (2.18) 

The terms of this double sum in which ai = a2 are written as a single contour 
integral and their sum is simply j^UA(ti+t2) of (p.ll|) . The Fourier transform 
of this term becomes 

-j^^ J J dhdhe-^'^'-^'^^UAih + h) = ^5(A - /i)p(A) (2.19) 
The remaining part, after the subtraction of the disconnected part, becomes 

1 r dudv _±_±_it^^_it^, 1 



J (27r?)2 (^u - y - - y ^ liA-j 

X IT 1 7 — r 1 ^ r 2.20) 

where the contours are taken around u = and v = a^. If we include also 
the contour-integration around the pole, v = u — j^, this gives precisely the 
term [/^(ti + ^2) of ( p.l9| ), which contributes to the delta-function part. This 
coincidence had already been noticed for the Laguerre ensemble [0. We 
are now in position to study the various properties of this random matrix 
problem with an arbitrary source A. 



3 Large N limit of the density of states 

For arbitrary A, the density of state p(A) was first found, in the large N 
limit, by Pastur[jl2|. The result may be easily recovered by summing dia- 
grams. Indeed in the large N limit the leading diagrams are planar, and 
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the one-particle Green function is the sum of "rainbow" diagrams. It fol- 
lows immediately that the self-energy is proportional to the Green function 
itself in the large N limit 0, and this leads at once to Pastur's result. From 
the contour-integral representation ( p.l2| ), let us show first how to recover 
this result . The average resolvent G{z) is written in terms of the evolution 
operator as 

dte-'''UA{t) (3.1) 
We substitute ( |2.12| ) for ^/^(t) and replace the product 

by its leading term in the large N limit, namely 

it ^ 1 

exp(-^ E ^) (3-3) 

If we define the density of states of the external matrix A 

1 ^ 

Po(a) = (3-4) 

a=l 

we may write this expression as 

exp(-2t / da^^) (3.5) 
J u — a 

Note that the "unperturbed" resolvent 
is related to po by 

Go{z) = [ da^ (3.7) 
J u + a 
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since Hq = —A. We obtain then easily 



' (3.8) 



dz J 27ri u + Go{u) — z 

We have now to specify the contour of integration in the complex u-plane. 
It surrounds all the eigenvalues of Hq and we have to determine the location 
of the zeroes of the denominator with respect to this contour. Let us return 
to the discrete form for the equation 

u + Go{u) = z (3.9) 

i.e. 

1 ^ 1 

u + -y: — = z (3.10) 

which posseses (A^ + 1) real or complex roots in the u-plane. For z real and 
large, N of these roots are close to the ej and one, which will be denoted u{z), 
goes to infinity with z as 

u{z) = z--^+0[^^) (3.11) 



Therefore, for large z, the contour encloses all the roots of ( p. 101 ) except 



u{z). When z decreases the contour should not be crossed by any other root 
of the equation, therefore it is defined by the requirement that only one root 
remains at its exterior. Therefore it is easier to calculate the integral (|3.8| ) 
by taking the residues of the singularities outside of the contour, rather than 
the N poles enclosed by this contour. There are two of them outside; one 
is u{z) and the other one is at infinity (since for large u, Gq{u) vanishes). 
Taking these two singularities we obtain 



dG _ ^ 1 

'dz ~ ~ 1 + #\ 

du{z) 

du{z) 



(3.12) 



The integration gives 



dz 

G{z)=z-u{z) (3.13) 
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(there is no integration constant since G{z) vanishes for z large; note that it 
does behave as it should as ^ for z large). This combined with (|3.9|) gives 
Pastur's self-consistent relation 

G{z) = Go{z - G{z)) (3.14) 



4 Universal correlations 

In the integral representation ( |2.20| ) we may neglect the terms t'^/N in the 
large N limit and replace the products as in (|3.5| ). This gives the large N-limit 

of f/jf^(ti,t2) as 



Noting that 



Hu{z,) - u{z,)] = -T^T-^-T-T^,-^-^ (4.2) 



dzidz2 ("^(-^i) ~ u{z2)y dzi dz2 

we obtain,through identical steps, the connected two-particle Green function 

^ liilti(.-i) - (4-3) 



This result was derived earlier by diagrammatic methods and was used 
to show that the singularity of the correlations, obtained when zi and Z2 
approach the real axis with opposite imaginary parts, is universal. 

However if we want to study the correlation function in the short-distance 
limit, we cannot use the resolvent any more (since we need to let the imagi- 
nary parts of ^1, Z2 go to zero before N goes to infinity). 

Returning then to (|2.2CI| ), and making the shifts, ti — »■ ti — iuN , and 
t2 ^ t2 — ivN, the two-level correlation function is remarkably factorized 



7=1 ^ "^7 + "/V 
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= -Kn{\iA2)Km{\2Ai) (4.4) 
This kernel K^i^ii A2) is further simphfied by the shift ti ^ t + if A^, 

J In J Imit ^ , i\(v — a.y) 



7=1 



Note that K]\f{Xi, Ai) reduces to the density of states. We replace again the 
product in ( |4.4| ) by its large N-limit, neglect ^ and integrate over t, leading 
to 

^ = -lm<f ^ - e--y (4.6) 

dXi 71 J 27rz u + Gq{u) — Xi + ie 

with y = N{Xi — A2). Therefore 

9Ai vr c/Ai 

= -^Ainife-?'"(^i-^^)) (4.7) 
Try c/Ai ^ ^ 

Since, from ( |3.13| ), 

u{Xi - ze) = Ai - ReG(Ai) - inpiXi) (4.8) 

we obtain 

K^iXi, A2) = -— e'^[^-^^^(^^)]sin[7ri/p(Ai)] (4.9) 
ny 

Repeating this calculation for K]\f{X2, Ai) we end up, in the large N, finite y 
limit, with 



Pc(Ai, A2) = -^sin2[7rp(^i±^)y] (4.10) 



Note that this result is independent of Ho (apart from the scale factor present 
in the densitty of states). In the case in which Hq vanishes it is also indepen- 
dent of the probability distribution of V|0]. It is natural to conjecture that 
this remains true for Hq non-zero as well, but we do not know of any method 
to prove it. 
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5 Conclusion 

The exact expressions for finite N of the correlation functions in terms of 

contour integrals, have allowed us to study the short-distance correlations 
for an arbitrary unperturbed Hamiltonian. We are not aware of any other 
method which could be used for solving this problem. The result is, as 
expected, a universal short distance behaviour, which depends on Hq only 
through scale factors. 
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